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Abstract 

In quantum systems with short-range interactions, causahty imposes nontrivial constraints on 
low-energy scattering parameters. We investigate these causahty constraints for systems where a 
long-range Coulomb potential is present in addition to a short-range interaction. The main result is 
an upper bound for the Coulomb-modified effective range parameter. We discuss the implications 
of this bound to the effective field theory for nuclear halo systems. In particular, we consider 
several examples of proton-nucleus and nucleus-nucleus scattering. For the bound-state regime, 
we find relations for the asymptotic normalization coefficients (ANCs) of nuclear halo states. As 
an application of these relations, we extract the ANCs of the excited 2"*" and 1~ states in ^^O from 
a— ^^C scattering data. 
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I. INTRODUCTION 



The constraints of causality for two-body scattering with finite-range interactions were 
first derived by Wigner pTj. The causahty bound can be understood as a lower bound on 
the time delay between the incoming and outgoing wave packets, At. When At is negative, 
the outgoing wave packet departs earlier than for the non-interacting system. However, 
the incoming wave must first reach the interacting region before the outgoing wave can 
leave. In low-energy scattering this manifests itself as an upper bound on the effective 
range parameter. In Ref. [2], Phillips and Cohen derived this bound for S-wave scattering 
with finite-range interactions. Some constraints on nucleon-nucleon scattering and the chiral 
two-pion exchange potential were considered in Ref. [3] , and relations between the scattering 
length and effective range have been explored for one-boson exchange potentials [4j and van 
der Waals potentials [5] . In Refs. [SI El the causality bounds for finite-range interactions were 
extended to an arbitrary number of spacetime dimensions and arbitrary angular momentum. 
The extension to systems with partial wave mixing was first studied in Ref. [8] . 

In this paper, we consider the causality constraints for the scattering of two charged 
particles with an arbitrary finite-range interaction. Our analysis is the first study of causality 
bounds that takes into account the long-range Coulomb force. The results presented here 
are relevant to studies of low-energy scattering of nuclei and nucleons using effective field 
theory (EFT). In particular, they are important for the application of effective field theory 
to nuclear halo systems (cf. Refs. P fTDj). 

This so-called halo EFT utilizes the separation of energy scales between the internal 
excitations of the core nucleus and soft halo physics. An effective Lagrangian is built order 
by order from local interactions and used to describe low-energy phenomena such as shallow 
bound states, charge radii, near-threshold resonances, radiative capture reactions, and soft 
photodissociation. Halo EFT has been used to describe neutron-alpha scattering [HI [12] 
and alpha-alpha scattering [13]. Bound single-neutron halo systems such as ^^Be [H] and 
^Li [151 [IS] have been studied as well as various two- neutron halo systems [T714T9] . 

There is an important connection between causality bounds and the convergence of ef- 
fective field theory calculations with increasing order [S] . For local contact interactions, the 
range of the interactions are controlled by the momentum cutoff scale of the effective theory. 
The cutoff scale plays an important role because the effective theory is typically not renor- 
malizable and new physics intervenes at higher energies. For calculations using dimensional 
regularization, the renormalization scale plays a similar role in regulating ultraviolet physics. 

The term "new physics", in this context, refers to details left out (integrated out) in the 
effective theory. In the case of halo EFT, these details are the finite size of the core nucleus 
and its internal excitations as well as the finite range of the interaction from pion exchange. 
Problems with convergence of the effective theory can occur if the cutoff scale is set higher 
than the scale of the new physics. However, it is useful to have a more quantitative measure 
of when problems may appear, and this is where the causality bound provides a useful 
diagnostic tool. For each scattering channel we use the physical scattering parameters to 
compute a quantity called the causal range, Rc- It is the minimum range for finite-range 
interactions consistent with the requirements of causality and unitarity. For any fixed cutoff 
scale, the causality bound marks a branch cut of the effective theory when viewed as a 
function of physical scattering parameters [8]. The coupling constants of the effective theory 
become complex when scattering parameters violating the causality bound are enforced. 
These branch cuts do not appear in perturbation theory; however, a nearby branch point 
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can spoil the absolute convergence of the perturbative expansion. 

Our results can be viewed as a guide for improving the convergence of halo EFT calcula- 
tions. In particular, if the cutoff momentum used in a calculation is too high, then problems 
with convergence may appear in some observables. Consequently, the causal range can be 
used to estimate the "natural" ultraviolet cutoff A of the effective theory as R~^- The nat- 
ural cutoff is optimal in the sense that no known infrared physics is left out of the theory 
and that all corrections involving the ultraviolet cutoff scale naturally [201 |2T] . Increasing 
the cutoff beyond the natural value will not improve the accuracy of the calculation. 

The causality bounds also have an impact in the regime of bound states. For two-body 
halo states, or more generally whenever there is a shallow two-body bound state close to 
threshold, the same integral identity that yields the causality bound for the effective range 
can be used to derive a relation between the asymptotic normalization constant (ANC) 
of the bound-state wave function, the binding momentum, and the effective range for the 
scattering of the two halo constituents. This relation can be shown to be equivalent to a 
result previously derived by Sparenberg et al. [22]. Its significance lies in the fact that the 
ANC is an important input parameter for the calculation of near-threshold radiative capture 
and photodissociation reactions. The causality bounds also constrain the range of model 
potentials that are fitted to scattering data in order to extract ANCs. 

The organization of this paper is as follows. We first review the general theory of scatter- 
ing for two charged particles with additional short-range interactions. Our analysis includes 
both attractive and repulsive Coulomb cases. In the next section we derive the charged- 
particle causality bounds for arbitrary values of the orbital angular momentum. Using 
the causality bounds, we extract and discuss the causal range for several nuclear scattering 
processes including proton-proton, proton-deuteron, proton-^He, proton-alpha, and alpha- 
alpha scattering. We then elucidate the relation for the ANC and extract the ANCs of the 
excited 2"*" and 1~ states in ^^O from a— ^^C scattering data as an application. We conclude 
with a summary of the main results and provide an outlook. In the appendices, we provide 
technical details of the derivation as well as numerical examples. 



II. TWO CHARGED PARTICLES WITH SHORT-RANGE INTERACTIONS 

In this section, we review some preliminaries that we will need later in order to derive the 
desired causality bounds. We consider a two-particle system with reduced mass /x interacting 
via a finite-range potential with range R. We write the interaction as a real symmetric 
operator with kernel F(r, r') satisfying the finite-range condition, 

V^(r,r') = \i r>R oi r'>R. (1) 

In the absence of Coulomb interactions the system is described by the radial Schrodinger 
equation, 

p'^Ui{r) = --T-^Ui{r) H — Ui{r) + 2/i / dr V{r, r') Ui{r') , (2) 

where by u^f^ we denote a solution for center-of-mass momentum p. Following the convention 
of Ref . [7] , the normalization of u^f^ is chosen such that for r > i? we have 

(r) = / [cot 6e{p) S,{pr) + C,{pr)] , (3) 
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where Sp and Ct are the Riccati-Bessel functions. 



If the particles carry electromagnetic charges Zic and ^26, respectively, there is a 
Coulomb potential in addition to the finite-range interaction. The radial Schrodinger equa- 
tion now reads 



P we{r) 



dr 



R 



d^ I" 7 

wt{r) + wi{r) + 2fi dr' V{r, r') We{r') + -Wi{r) (4) 



with the Coulomb parameter 



The normalization of w^f^ is now chosen such that for r > R we have 



7 = 2/i ■ aZiZ2 . 



(5) 



w'f^ (r) = p'Cr,,e cot 6e{p) Ff> (r) + G'[> (r) 



(p)/ 



(6) 



where P'f'' and G^^^ are the regular and irregular Coulomb wave functions [23] , respectively, 
and 81 is the phase shift of the full solution wf^ compared to the regular Coulomb function 
P'f^ [21]. 1 The factor C^,^ is given by 



7] 



C. 



1_ 
27rrj 



and 



22e 



[{2i + 



(7) 
(8) 

(9) 



A. Coulomb wave functions 

For the Coulomb wave functions F^-^^ and Cf'^ we use the conventions introduced by Yost, 
Breit and Wheeler in Ref. [23]. Explicitly, we write them as [251 EE] 

Gf\r) = |^^f^^~i| e-'f(i-+-)iy..^(^) + iFi'\r) , (10b) 
\\. + m — K,)\ 

where ^ 

p = p ■ r , z = 2ip , n = 17] , m = ^H — . (11) 

^Note that for f = our normalizatfon is the same as chosen m Ref. |21], ie., for r > i? our solutfon w'"^'^ 
coincides with the function if defined in Eq. (42) of that paper. 



r(2m + 1) 



(10a) 
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The functions M^^m and WK,m are Whittaker functions, which can be expressed in terms of 
hypergeometric functions as 

M^^miz) = e-5^z5+-,Fi (i + m - 1 + 2m; 2) , (12) 

W^^Uz) = e-5^z5+-f/(l + rn- K,l + 2m;z) . (13) 
iFi(a, b; z) is Kummer's function of the first kind, 

iFi(a,6;^) = ^^-y^ , a^") = a(a + 1) ■ ■ ■ (a + n - 1) , (14) 



n=0 



and 



^(«, ^) = T^7?I— ^) + - 6 + 1, 2 - 6; z) . (15) 

1 (1 + a — 0) 1 (a) 

A comprehensive discussion of the functions F^\r) and G^/\r) can be found in Ref. |27] . 
For asymptotically large p = pr, the Coulomb wave functions behave as 

) (r) ~ sin(p - £7r/2 - r/ log(2p) + (x,) , (16) 
(r) ~ cos(p - in/2 - rj log(2p) + a^) (17) 

with the Coulomb phase shift 

ae = aigTii + 1 + 17]) . (18) 
In the limit p — )■ 0, their behavior is 

F,"''(r)~C,,,p'+'. (19a) 



B. EfTective range expansions 

For a system without Coulomb force and only a short-range interaction, one has the 
well-known effective range expansion 

p^'-^' cot 5eip) = -- + lr,p'' + -- - , (20) 
ae Z 

where a£ and are the scattering and effective range parameters, respectively. For a system 
of charged particles we have instead the Coulomb-modified effective range expansion 

Cl,p''^' cot Up) + 7 hip) = -i, + irf / + • • • , (21) 

where 

heip) = P''%^ h{r]) , (22) 

h{ri) = Reipiirj) — log |?7| , (23) 
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and V'(^) = r'(2;)/r(2;) is the logarithmic derivative of the Gamma function, also called 
digamma function. For i = 0, the expansion simplifies to 



CloP cot 6o{p) +-fh{r]) 



'-roP^ + 



(24) 



A derivation of Eq. (24) for the case of proton-proton scattering can be found, for example, 
in Ref. [21]. See also Ref. [28] for a detailed discussion. The analytic properties of the i = 
modified effective range function are investigated in Ref. [29J.^ 

In Ref. [30], Bolle and Gesztesy derived a very general form of the Coulomb-modified 



effective range expansion for an arbitrary number of spatial dimensions. Specializing their 



result to the three-dimensional case, a version of Eq. (21) can be written as 

1 



Cl,p''^UcotUp)-^)+l~hi{p) 



1 



(25) 



wi 



ith^ 



{2pf' 




r(2£ + 2)2 


r(l + ir7)|2 



The latter function can be rewritten using 



^2 



221 


r(£ + l + ir7)|2 


r(2£ + 2)2 


r(l + ir/)|2 



IVf] 



(26) 



(27) 



with Cl Q as defined in Eq. ([s]). The expressions given here reproduce Eqs. ( |2T| ) and ^ when 
one explicitly assumes that the momentum p is real. In fact, one has to rewrite Eq. (26) in 



this manner in order to get an effective range function that is analytic in around threshold 



III. CAUSALITY BOUNDS FOR CHARGED PARTICLES 

With the Coulomb wave functions and Coulomb-modified effective range expansions at 
our hands, we can now closely follow the derivation presented in Ref. [7] for scattering in 
the absence of Coulomb interactions. 



^As a remark we note that on first sight the expansion given in Eq. (51) of Bethe's paper [21] seems to be 



different from the one given here in Eq. ( 24 1 , which is the same as given in later publications referring to 
Bethe's result. The ?7-dependent function on the left hand side of Bethe's expansion appears to differ from 
our h{ri) by two times the Euler-Mascheroni constant ^e- This apparent conflict can be resolved by noting 
that the 5(77) in Eq. (51) of Ref. (SJj is not the function defined in Eq. (47a) of the same paper, but rather 
given by lim^j_j.oo [3(77) — gijji)], where in this latter expression the g from Eq. (47a) is meant. The limiting 
process then yields exactly the term —2^e- 

■^This definition essentially comes from combining Eqs. (4.1) and (4.2) in Ref. jSU], with the correction that 
the exponent in Eq. (4.2) should be —2 rather than 2. 



6 



A. Wronskian identities 

We consider solutions of the radial Schrodinger equation Q for two different momenta 
Pa and pb- Introducing the short-hand notation 

WA,B{r)=wt^^''\r), (28) 
i.e., suppressing the angular-momentum subscript £ for convenience, we get 



(p| - Pa) I dr'wA{r')wB{r') = {wbw'a - waw'b)\1 



R 



— 2/i / dr' / dr [wB{r)V{r,r')wA{r') — WA{r)V{r,r')wB{r')] (29) 



by subtracting wa times the equation for wb from that for wb multiplied by wa, as it is 
done in Ref. [7], and integrating from a small radius e to r. 

We assume that our interaction V{r,r') is such that it alone permits a solution that is 
sufficiently regular at the origin, i.e., 'U^(O) = and drUe stays finite as r — i- 0, where 
is a solution of Eq. As boundary condition for the solutions wa,b of the full radial 
Schrodinger equation we can then demand as well that they vanish with finite derivative 
at the origin. If we only had the Coulomb potential and no additional interaction, this is 
fulfilled by the regular Coulomb function F^\r), of. Eq. (19a). We can thus take the limit 
e — > in Eq. ( 29 ) and get the relation 

pr 

W[wB,WA]ir) = (pI-pD / dr' WAir')wBir') , (30) 

where the Wronskian is defined as 

W[wB,WA]{r) = WB{r)w'A{r) - WA{r)w'B{r) . (31) 



B. Rewriting the wave functions 

Following further the derivation presented in Ref. [Tj, we re-express the solutions w^^\r) 
in terms of functions f{p, r) and g{p, r) such that 

w^;\r) = p^'^'Cl, cot Up) fiP, r) + 9{p, r) (32) 
for r > R, with f{p,r) analytic in p^, 

fip,r) = foir) + Mr)p' + Oip'), (33) 

and 

g{p,r) = g{p,r) + (j){p) ■ f{p,r) . (34a) 

The g{p,r) contains a term which is non-analytic in p"^ and is proportional to f{p,r). The 
remainder g{p,r), however, is analytic in p"^, 

~g{p,r) = go{r) + g2{r)p' + 0{p^) . (34b) 
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Explicit expressions for these functions can be obtained from the results of Bolle and 
Gesztesy [30]; they are discussed further in Appendix [Bj Combining Eqs. ^ and (32), 
we find 

fiP,r) = ^^F^'\r) (35a) 

and 

g{p,r)=p'C,,eGf\r). (35b) 



We insert now the modified effective range expansion (21) into the asymptotic Coulomb 



wave function ( 32 ) . Note also that 



(36) 

where the last step follows from the relations in Appendix |Bj The term involving h£{p) 
exactly drops out and we are left with 



1 1 



— ^ + ^rfp^H ] f{p,r) + g{p,r) , r>R. 



(37) 



Thus it is possible to choose a normalization such that w^f\r) is analytic in p^. Combining 
this with the expansions (33) and (34), we arrive at 



wPir) = — ^/o(r) + go{r) + p^ 



^ do 



+ 0{p') . (38) 



C. Causality bounds 

From here, we can proceed exactly as in Ref . [7] . For the Wronskian of two solutions wa 
and wb for r > R we find 



1 



W[wB,WA]ir) = (p^ -p1)<; -r'^W[fo,go]ir) 



W[f2,fo]{r) 



- -c mf2,go]{r) - W[g2, /o](r)] + W[g2,go]{r) } + O^pX^) . (39) 



Note that in the 0{p\q) we have also included terms of the form p^Ps- We set p^i = in 



Eq. (30) and furthermore take the limit pb 0. Using the expansion (39), we get 

2 



/ dr' 




Jo 





(40) 



for r > R, with wf^\r) = \im.p^owf' (r) and 



ip)i 



6f (r) = 2W[g2,goW) - ^ {W[f2,go\{r) + W[g2, /o](r)} + -^W[h, /o](r) . (41) 
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Written as a function oi p = p ■ r, the Wronskian of the Coulomb wave functions is^ 

W[Fi'\ G?]{p) = -W[Gf\ = -1 . (42) 

Since d/dr = p ■ d/dp and W[f, f] = 0, we also have 

W[f,g]{r) = W[f,~g]ir) = -l. (43) 



Plugging in the expansions p3| ) and (34), we see that iy[/o, fi'o]('^) = —1 for the leading-order 
functions, and W[f2, go]{r) = W[g2, fo]{r) for the terms at O(p^). Inserting these relations 
into Eq. (40), we get 

r? = b?(r) - 2 



/ dr' 


wl '{r) 


Jo 





where bf{r) has been simplified to 



bfir) = 2W[g2,go]{r) 



-,c 



W[h,go\{T) 



C^2 



Since the integral in Eq. (44) is positive definite, the resulting causality bound is 

rf < &f (r) , V r > i?. 



(44) 



(45) 



(46) 



D. Calculating the Wronskians 



We now want to derive the explicit form of the function h'^{r). To do this, we need 
expressions for the Wronskians that appear in Eq. (45). We can obtain them by first noting 



that f{p,r) and g{p,r), being linear combinations of Coulomb wave functions (with p- 
dependent coefficients), are solutions of the Coulomb Schrodinger equation. 



dr" 



x{p, r) = 



(47) 



which, of course, corresponds to setting V{r,r') = in Eq. (|4]). Here and in the following, 
X stands for either / or g. Inserting the expansion 



x{p, r) = Xo{r) + p^ X2{r) + O(p^) 
into Eq. (47) and comparing orders in p^, we find that 

d2 



dr" 



+ 



+ 1) , 7 



0. 



(48) 



(49) 



i.e., Xq is a solution of the zero-energy Coulomb Schrodinger equation, and 

d2 



dr" 



+ 



+ 1) , 7 



X2(r) = Xo(r) . 



(50) 



*See, for example, Eq. (14.2.4) in Ref. [3T]- 
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From this we readily obtain the differential equations 

^W[h, fo]ir) = [foir)]\ (51a) 

±W[g2,go]{r) = [go{r)]\ (51b) 

^W[f2,go]ir) = foir)go{r) (51c) 

for the desired Wronskians. Put together, this yields a simple first-order differential equation 
for bf{r), 

^&.'^W = 2(^7o(r)-^/o(r)y. (52) 

From Eqs. (A. 7) and (A. 8) in Ref. PDJ we have the explicit expressions 

(21 + IV 

/o(r) = v^/2,+i(2y^) , (53a) 
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21+1 



9o{r) = 7^2^+1 (2 v/^) (53b) 



for 7 > 0, where la and Ka are modified Bessel functions, and 



/o(r) = -^^^ J,e+i{2V^) , (54a) 
9o{r) = ^^/^ Ij, v^l^2.+i(2v^) (54b) 



for 7 < 0, where Ja and Y^, are the ordinary Bessel functions. ^'^ 



Using these expressions for /q and go and Eq. (52) we can determine bf{r) up to an inte 



gration constant. In order to fix this constant, we must work directly with the Wronskians 



in Eq. (41). Before we do that, however, we first discuss the general form of bf{r). Let us 
break apart the function as a sum of two functions, X£{r) and Yi{r), and a constant term 

}P,{r)=X,{r)^Y,{r)^Z,. (55) 

We take X({r) to be a function consisting entirely of a sum of terms that have a pole at 
r = 0, ranging from order 1 to 

X,[r) = Y,^Lmr-"^ ■ (56) 

■m=l 



^From Eq. (9.1.50) and the remark above Eq. (9.6.41) in Ref. [3T] it is clear that these /g and are indeed 



solutions of (49 1. 

^BoUe and Gcsztesy actually give an expression for in the attractive case (7 < 0) that involves the Hankel 
function ij'^^ times i instead of Y . With that, however, gg would not be real, which it should be. Our gg as 



in Eq. ( 54b I is taken from the results of Lambert 
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By furthermore requiring the function Yi{r) to vanish at r = 0, the decomposition in Eq. (55) 
is unique. 

Where exactly the contributions to the three terms in the decomposition come from can 



be inferred from the behavior of f{p,r) and g{p,r) at the origin. From Eq. (19a) combined 
with Eq. (35) we find that 

f{p,r) ~ as r 0. (57) 
This imphes that every term in the expansion of f{p,r) is (9(r^"*"^). Therefore, 



hmiy[/2,/o](r) = 

r— ^0 



(58) 



for all i, which means that this Wronskian only yields contributions to the Y£{r). 

Furthermore, from Eq. (17) in Ref. [23] we know that the irregular Coulomb wave function 
has the asymptotic behavior 



G 



(p) 



D^ri,e P ^ as p -T- - 



(59) 



with D^^n fulfilling Cn^iDn/ = 2£ + 1. Using Eq. (35b) then yields 



g{p,r) 



2i + l 



as r — )■ . 



(60) 



We note that go{r) has exactly the same behavior near r = and can thus show that g2{r) 



is subleading as r — 0, g2ir) 



-£+c 



for c > 0. From this we infer that 



hmW[fo,g2]{r) = 

r— >0 



(61) 



for all i, so also from this Wronskian we only get contributions to Y£{r). Both the singular 
X£{r) and the constant Ze, therefore, only come from the Wronskian W[g2, go]{r). 
For i = the situation is still simple because the above analysis also tells us that 



\imW[g2, go]{r) =0 for i = 0. 



r-s>0 



(62) 



i.e., 6f (r) is given entirely by Yo{r). With the knowledge that it vanishes at the origin, it is 
actually straightforward to give an explicit expression for b^{r) in terms of antiderivatives 
of the right hand side of Eq. (52), where one has to insert the /o('") and go{r) from Eqs. (53) 
and (54). The result, obtained by integrating from to r, is 



2r^ -9 /3 



4^2 



K)"'G?;U47r 



1 

0,1,-2 



+ ^^-ir'^ ( 47r 



-1 



' 2 



-1,0,1,-2 



for the repulsive case, and 



(63) 



2^3 



1F2 ( ^; 2, 4; 47r ) + ^^r-" (a^) ' Gf^ ( -47r 



1 _i 
2 ' 2 
0,1,-2,-1 



+ 27r^ 



^iF2('|2,4;47r 



27r' ^3,1 



G 



3,5 



-47r 



-1 i 
' 2 ' 2 

-1,0,1,-2,- 



(64) 
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for an attractive Coulomb interaction. In the above equations, pFq and denote the 

(generahzed) hypergeometric and Meijer G-functions, respectively. 

For general i > 1, W[g2, go]{r) is singular at r = and the analysis becomes more 
complicated. For practical purposes, one can simply use power-series expansions for the 
Bessel functions that appear in the expressions for the zero-energy functions and integrate 
these term by term until a desired precision is reached. The only additional ingredients 
needed are the values for the constant terms Ze because these are obviously not generated 
by the integration. In Table |T] we list these constants for i = 0,1, 2. 



i 





1 


2 







7 




^3 


{ 79 7£\ 


^6 9 ; 


*v 21600 360 y 



TABLE I: Constant term in Eq. (55) for ^ = 0, 1, 2. = 0.577216 ... is the Euler-Mascheroni 



constant. The values are the same for repulsive (7 > 0) and attractive (7 < 0) Coulomb potentials. 



In Appendix |A] we describe how the values in Table |T] can be obtained with the help 
of computer algebra software. In that appendix we also give explicit expressions for the 
causality bound functions 6f (r) and 6^(r). 



IV. EXAMPLES AND RESULTS 
A. The causal range 



The causality bound (46) can be rewritten as 

6f(r) -rf > Vr > i?. (65) 
In cases where the details of the interaction (in particular its range, assuming that a de- 



scription with finite-range potentials is applicable) is not known, one can use Eq. (65) to 
define the causal range of a scattering system as that value of r for which the bound is 
just satisfied, i.e., 

bf{R,)-rf = 0. (66) 



We note from Eq. (52) that the derivative of bf{r) is non- negative, 

^bf{r) = 2(^go{r)-^fo{r)y >0. (67) 

Hence, bf is an increasing function of r and the causal range is defined uniquely. For the 
case that Eq. (66) does not have a solution {i.e., if bf{r) is positive already for r = 0), 
we define the causal range to be zero. The importance of the causal range is given by the 
fact that it can be interpreted as the minimum range a potential is allowed to have to be 
consistent with causality. If for a given system the values in individual partial waves differ 
significantly, the maximum value should be taken as the causal range of the underlying 
potential. Alternatively, one can model the interaction with an ^-dependent potential. For 
effective field theories with short-range interactions such as halo EFT, the causal range 
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constrains the allowed values of the momentum-space cutoff or the lattice spacing used in 
numerical calculations. 



In Fig. [T]we plot the left-hand side of Eq. (65) as a function of r for the case of proton- 
proton S-wave scattering. The causal range can then be read off as the point where the 
function becomes zero. Fig. [2] shows analogous plots for a-a S- and D-wave scattering. In 
this system, there are visible error bands due to the larger uncertainties in the effective range 
parameters. 




FIG. 1: Causal range plot for S-wave proton-proton scattering. 
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FIG. 2: Causal range plot for S-wave a-a. scattering. 

In Table |TT] we give a summary of the causal ranges that one finds for various two- 
body systems of light nuclei where low-energy scattering parameters and/or phase shifts are 
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available from experiments. The results are briefly discussed in the following subsections. 



System 


Reference 


Channel 


a^/ fm^^+i 




Causal range / fm 


p-p 


[33] 




-7.828 ±0.008 


2.80 ±0.02 


1.38 ±0.01 


p-p 


[31] 




-3.03 ±0.11 


4.22 ±0.11 


2.33 ±0.05 


p-p 


[34J 


'Pi 


2.013 ±0.053 


-7.92 ±0.17 


« 0.03 


p-d 


[35J 


2^1/2 


2.73 ±0.10 


2.27 ±0.12 


3.90 ±0.15 


p-d 


I3H] 


'^'^1/2 


4 


-2.8 





p~d 


[35] 


4 e 

'-'3/2 


11.88 ±0.40 


2.63 ±0.02 


2 20+°-°'^ 
— 0.06 


p~d 


[36J 


4 e 

'-'3/2 


11.11 


2.64 


2.29 




i37j 


5*0 


11.1 ±0.4 


1.58 ±0.12 


1 qo+0.21 




m 




9.04 ±0.14 


1.50 ±0.06 


1 97+0.10 

-•^•^'-0.09 


p-a 


m 


-51/2 


4.97 ±0.12 


1.295 ±0.082 


1 09+0.40 

±.OZ_Q 21 


p-a 


[38] 


Pl/2 


-19.36 ±0.50 


0.349 ± 0.021 


2.65 ±0.07 


p-a 


[38] 


P3/2 


-44.83 ±0.51 


-0.365 ±0.113 


'-'•^^-0.10 


a-a 


[13J 


s 


(-1.65 ±0.17) • 10^ 


1.084 ±0.011 


9 f^8+0.14 
^•^°-0.13 


a— a 


[39J 


D 


(-7.23 ±0.61) • 10^ 


(-1.31 ±0.22) • 10^3 


9 nQ+0.09 
^•^■-•-0.07 



TABLE II: Summary of causal-range results obtained from experimental input for various few- 
nucleon systems. 



1. Proton-proton scattering 

For p-p S-wave scattering one flnds a causal range of about 1.38 fm. This value is very 
close to the value obtained by assuming that the typical length scale of the N-N interaction 
is set by the inverse pion mass, Hc/Mt^ 1.4 fm. The value one flnds in the '^Po-channel is 
somewhat larger [R^ ~ 2.3 fm), whereas the ^Pi effective range parameters impose almost 
no constraint on the range of the nuclear potential in this channel. This suggests some 
signiflcant differences in the radial dependence of the interactions for the ^Pi channel. In 
fact, the tail from the one-pion exchange interaction is repulsive in the ^Pi-channel while it 
is attractive in the ^Po-channel (c/. the discussion in Ref. [8]). 

For effective fleld theory calculations with purely local interactions {i.e., pionless effective 
fleld theory), our results suggest to keep the cutoff momentum smaller than M^- for the ^Sq 
and ^Pq channel. However there is more freedom to take a higher cutoff for the ^Pi channel. 

In Ref. [8] , causality bounds were investigated for neutron-proton scattering. The results 
(Pc = 1-27 fm for ^6*0, Rc = 3.07 fm for ^Pq, and R^ = 0.23 fm for 'Pi) are qualitatively very 
similar to what we flnd for the p-p system. This would indicate only a moderate amount of 
isospin breaking. 

2. Proton-deuteron scattering 

There are several experimental determinations of p-d effective range parameters. In 
Table [lH we have included results from Arvieux [35j and Huttel et al. |36j- While for the 
quartet-channel there is a good agreement between the scattering lengths and effective ranges 
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(and, of course, of the resulting causal ranges, which come out as 2.2 — 2.3 fm), there is a 
large discrepancy for the doublet-channel results. 

The difficulty of determining the proton-deuteron doublet-channel scattering length has 
previously been discussed by Orlov and Orevkov |1D]. Comparing different models, the 
authors conclude that a*-^ = 0.024 fm is currently the best theoretical estimate for the 
doublet-channel p-d scattering length. From Table 3 in Ref. |10] one reads off that the 
corresponding value for the Coulomb-modified effective range is as huge as r*^ = 8.23-10^ fm. 
Inserting these numbers into the causal-range calculation one gets a very large value of 
Rc ~ 8.15 fm. As a consequence, no definite conclusion can be reached with the currently 
available data. 



3. Proton-helion scattering 

For the scattering of protons off a helium nucleus we were able to find data for both p-^He 
and p-a scattering. In the first case, there was only enough data available to calculate the 
causal range for the S-wave channels. Since both the scattering lengths and effective ranges 
are very similar for the singlet and the triplet channel, so are the resulting causal ranges, 
which come out as approximately 1.3 fm. 

Incidentally, one finds almost the same value for the S-wave in p-a scattering. For this 
system, it is interesting to compare to the neutron-alpha system, where there is no Coulomb 
repulsion in the scattering process. Results for the n-a causal ranges can be read off from 
Fig. 5 in Ref. [7j (obtained using effective range parameters from Ref. [3H|). Even though 
from the plot one only gets quite rough estimates, one clearly sees that the results for the S1/2 
and Pi/2 channels agree very well between p-a and n-a scattering, which as in the nucleon- 
nucleon case discussed above could be interpreted as only a moderate amount of isospin 
breaking. However, the causal ranges for the P3/2 channels are very different (~ 0.5 fm for 
p-a, ~ 2 fm for n-a). It is an interesting question if this discrepancy hints at an error in 
the extraction of the effective range parameters (either for one of the systems or possibly 
both), or if there actually is a physical reason behind the difference in the causal ranges. 



4- Alpha-alpha scattering 

For a-a S-wave scattering we use the values given by Higa et al. (see Ref. [13] and 
experimental references therein), = (—1.65 ±0.17) -10^ fm and Tq = (—1.084 ±0.011) fm 
to find a causal range of about 2.58 fm. 

For the i = 2 channel, no effective range parameters could be found in the literature. 
We have thus used the phase shift data collected in the review article by Afzal et al. [39] 



to perform the fit to the effective range expansion (Eq. (21) with i = 2) ourselves. By 
including the phase-shift data up to i?iab ~ 6.5 MeV we find = (—7.23 ± 0.61) ■ 10^ fm^ 
and = (—1.31 ± 0.22) ■ 10^'^ fm^'^. However, the fit is strongly dominated by the O(p^) 
shape parameter, so the actual uncertainties of and should probably be somewhat 
larger. For the causal range in this channel we find a value of about 2 fm, which is just 
slightly smaller than the S-wave result. 
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V. RELATION FOR ASYMPTOTIC NORMALIZATION CONSTANTS 



We now discuss the connection between the integral relations we have derived above and 
asymptotic normalization constants (ANC), denoted in the following as A. They are defined 
by writing the bound-state wave function in the asymptotic region as 

Weir) = A ■ Vr_i^_^+i(2Kr) for r>R, (68) 

where w^(r) is a bound-state solution (with momentum p = in, k > 0) of the Schrodinger 
equation Q, normalized according to 

X) 

dr\we{r)f = 1. (69) 



is the Whittaker function defined in Eq. (13). At very low-energies, the radiative 
capture cross section to a two-body halo nucleus is proportional to \A\'^. 



Derivation of the relation 



Our starting point is Eq. (44) derived in Sec. Ill, from which we already obtained the 
causality bounds for the effective range parameter. For bound-state momenta, we write it 

as 



IT 

I dr' 


w} '{r) 


Jo 





for all r > R. 



(70) 



with bf{r) as defined in Eq. (45). In order to derive from this a relation for the ANC, we 
need the precise connection between the radial wave functions ti'^"^^ and the solution Wi{r) 
used in Eq. (68). Formally, we can write 



cot6AiK)Fi'^\r) + G 



for r > R. 



(71) 



but we need to be careful with the way the analytic continuation to the bound-state regime 
is done. We start with the assumption that cot 6i{iK) = i for a bound state and thus consider 
the linear combination G'f\r) + iF^^\r), not explicitly setting p = ik for the moment. A 
short calculation shows that 



(r) + iF^\r) = e^^'e-'W+-^^) W_,^^^^{-2ipr) 
with the Coulomb phase shift defined via [321 SI] 

T(£ + l + ir7)' ' 



r(£+l-ir7) 



(72) 



(73) 



Note that ir] will be real for bound-state momenta. In the following, we only consider the 
case of a repulsive Coulomb potential (i.e., 7 > 0), so ir/ is in fact a positive number. For 
the Cr,,i, we use the expression |42] 



2^e-'^ [T{i+l + i7])T{i + l-i7])] 
r(2£ + 2) 



(74) 
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cf. also Eq. (|B2|) in Appendix |B[ Combining Eqs. (|73|) and (|74|) we find that the problematic 

(75) 

2^T{e + l + ir]) 



terms T{i +1 — 17]) and exp{'n-ri/2) drop out, and we can finally write 

wf\r) = K^C'^,,iy_i^^,+ i(2Kr) for r>R 

where we have defined 



r(2£ + 2) 



Comparing with Eq. (68) we readily infer that 



K 



Using the overall normalization of the wave function, we can rewrite Eq. ( 70 ) as 

2^ /"OO r , , -,2' 



br(r) - 21im 



K 



Ct,- 



dr' 



w) '[r 



, r > R. 



(76) 



(77) 



(78) 



Canceling the r -dependence 



At this point it is important to note that Eq. ( 78 ) can be rigorously valid only for zero- 
energy bound states, i.e., if one strictly considers the limit k — )■ 0. Since the left hand 
side is a constant, the r-dependence of bf{r) has to cancel that of the integral for any 
r > R. In the scattering regime this is ensured by taking zero-energy limit, i.e., considering 
scattering directly at threshold. In the bound-state regime, on the other hand, one is of 
course interested in finding a relation for the case where there is a bound state close to 
threshold, i.e., where k is small but finite. It is obvious that in principle this poses a 
problem because the bf{r) in Eq. (78) is a strictly increasing function of r, whereas the 
integral is always bounded and in fact becomes smaller when r is made larger. 

However, up to corrections of higher order in k, the cancellation still works, i.e., the 
leading r-dependence drops out. In Ref. |13], where we derived an analogous ANC relation 
for the case without any long-range force, we established the cancellation by carrying out 
the integral over the asymptotic wave function (a Riccati-Hankel function, in that case) 
analytically and then expanding the result in in powers of k. To show it here in a more 



general way, we recall from Eq. ( 38 ) that 



1 



fo{r)+go{r) + 0{p' 



(79) 



It is straightforward to assume that this is valid also in the bound-state regime since we are 
working with a wave function explicitly analytic in p^. Hence, 



fo{r)+go{r)] + 0{k 



2d^ 



bf{r) + 0{k' 



(80) 



where in the last step we have inserted Eq. (52). This directly tells us that b^{r) in Eq. (78) 



cancels with the integral up to higher-order terms and, possibly, an integration constant. 
This situation is already familiar from the calculation of the Wronskian W[g2,go\{T). Again 
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the constant term can only come from the integral over goir)"^- To determine it, we note 
that 

/•oo 

t2 



dr'[go{r')Y < oo if r > 0, (81) 
i.e., the integral is actually convergent and gives a Meijer G-function, as already encountered 



in Sec. Ill D The constant term coming from this function has to be combined with the 
already known Zi in bf{r). We define 



(82) 



/■oo 

AZe = W[g2,go]ir)+ / dr' [goir')] 

J r 

which is explicitly r-independent, to finally arrive at 



or, equivalent ly. 



\A\ 



a 



AZi + 0{k^ 



-1/2 



(83) 



(84) 



both valid in the limit where k — )■ 0. In Table |III| we give the resulting values for the AZ^ 
for £ = 0,1, 2. 








1 


2 




-1/(37) 


-7/108 


-I773/IO8OO 



TABLE III: Integration constant AZ^ in Eqs. (83 , 84) for £ = 0, 1, 2 



Finally, we note that Sparenberg et al. [22] have previously derived an ANC relation 
equivalent to ours, but written in terms of the scattering length instead of the effective 
range parameter. The equivalence of the two relations up to the contribution of higher 



order shape parameters is shown in Appendix D 



B. Application to the oxygen-16 system 

The ^^O nucleus has two excited states lying just below the a-^^C threshold, a 2+ at about 
—245 keV and an even more shallow 1^ at only —45 keV. The properties of these states 
play an important role in astrophysical helium burning processes [HHISI- In the following, 
we calculate asymptotic normalization constants for both states under the assumption that 
they can be described in a a + ^^C halo picture. 

We use the recent data obtained by Tischhauser et al. |17] (for the actual phase shifts see 
their Ref. [32]) in order to extract the Coulomb-modified effective range for the a-^^C P- and 
D-wave channels. Focusing first on the D-wave, we note that the combination of the strong 
Coulomb repulsion and the £ = 2 centrifugal barrier makes the low-energy phase shifts very 
small over a wide energy range. Moreover, there is a narrow resonance at a center-of-mass 
energy of about 2.7 MeV, which strongly constrains the energy region for a straightforward 
fit to the effective range function. From a simple fit to the data up to the narrow resonance 
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one only obtains an effective range parameter with an uncertainty too large (50%) for an 
extraction of the ANC since the latter depends on r'^ very sensitively. 

To mitigate this problem, we use the position of the 2"*" oxygen state as an additional 
input parameter for a self-consistent extraction of the effective range. In the following, we 
describe this procedure in more detail. 

At the position of the bound state, where cot6£{iK) = i, the Coulomb-modified effective 
range expansion ( 25 ) reduces to 



-r9 + 



(85) 



where now 



(2iK 



\2e 



(2£+l)!= 



n( 

s=l 



^(i^) + 7^- log(i^/) 

ZIT] 



A straightforward calculation shows that for the prefactor we have 



7 

2i«; 



(86) 
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(2£ + l)!' 



1 



2t 



(2£+l)!^ 



(87) 



Furthermore, for the digamma function we have the asymptotic expansion (Eq. (6.3.18) in 
Ref. [31]) 

1 °° R 



2z 



n=l 



for \z\ — )■ oo, where the are the Bernoulli numbers, 



5n 



1, B, 



1 

2' 



1 

6' 



(89) 



Noting that the sum in Eq. (88) only starts at n = 1, we find that 7/1^(1^) only starts at 
order k^, and inserting the precise relation into Eq. (85), we arrive at 



1 



2 ^ 3{2i 



1 



(90) 



Neglecting the effective range contribution at leading order, we insert the binding mo- 
mentum K ^ 0.187 fm^^ of the 2+ state to get a first approximation for the scattering length 
parameter. This is then used to constrain a subsequent fit to the phase shift data up to 
about 2.6 MeV, i.e., just below the narrow resonance. We include a single 0{j)'^) shape 
parameter in the fit since some curvature is clearly necessary. The obtained from this is 
then used in Eq. (90) to get a better result for 1/a^, which, in turn, is fed back into the fit. 
Iterating this procedure a couple of times yields a well- converged self-consistent result for 
r^. After eight iterations, we find 



(5.94 ± 0.35) • 10'^ fm" 



(91) 



for the a-^^C D-wave. Including a second shape parameter in the fitting procedure only 
changes this result within the given uncertainty, so we conclude that for the energy range we 
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have been fitting, a single shape parameter really is sufficient to account for the curvature. 
Inserting the fit result into Eq. (Islj) yields 



|v4(2+)| = (2.41 ± 0.38) ■ 10^ frn'^/^ (92) 



for the 2+ state in O. Including in Eq. (83) an 0{k ) term of the order of the shape 



parameter gives a consistent result for |74(2''")| within the error given above. 

When we apply the same same procedure to the 1~ state just 45 keV below the a-^^C 
threshold, we also see a nice convergence and obtain the results 

= (4.546 ± 0.002) • 10"^ fm"^ (93) 

and 

\A{l-)\ = (1.188 ± 0.024) ■ 10^^ fm-^/' . (94) 

The overall picture for the 1~ state is somewhat more complicated, however. On the one 
hand, the extraction of the effective range parameters is easier in this case because there is 
no narrow resonance limiting the fit range (we have used the data up to i?cm = 3.75 MeV). 
By doing a simple fit without the self-consistent iteration we get values for the effective range 
and the ANC that are slightly smaller, but overlap with the results given above when the 
respective uncertainties are taken into account. On the other hand, allowing for a second 
shape parameter changes the value of the effective range quite dramatically to 0.046 fm~^, 
which leads to an imaginary ANC. This could indicate that the effective range parameters 
from the simple fit violate the causality bound. The fact that has to be such that 



Eq. (83) yields a real value for the ANC can be interpreted as a weaker remnant of the 
original causality bound relation. Alternatively, the cluster picture might not be applicable 
for the shallow 1~ state in ^^O. 

Our result for the ANC of the 2"*" state is about a factor five smaller than the value 
\A{2+)\ = (1.11 ±0.11)- 10^ fm-^/2 obtained by Brune et al. [Mj , while the value for \A{l-)\ 
is about a factor two smaller. Other, more recent determinations [l8l HH] have found even 
larger values. Sparenberg et al. [50] have carried out a similar analysis for the 2+ state 
based on their ANC relation (c/. Ref. [22]). They subtracted the narrow resonance from the 
a-^^C D-wave phase shift data in order to extract a set of higher-order shape parameters 
and concluded that present day data are not sufficient to constrain the ANC strongly. As 
discussed above, our approach of performing a self-consistent ffi to the data below the 
resonance constrained by the separation energy and using the effective range instead of the 
scattering length as input in the ANC relation improves the stability of the extraction. 
However, compared to other determinations our ANC values are generally smaller. This 
issue requires further study. 



VI. CONCLUSION 

In this paper, we have investigated the constraints imposed by causality on the low-energy 
scattering parameters of charged particles interacting via a short-range interaction and a 
long-range Coulomb potential. Similar to the case of neutral particles without Coulomb 
interaction [2]-(8], our considerations lead to a constraint on the maximum value of the 
Coulomb-modified effective range. 
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While conceptionally straightforward, the calculation of the Wronskians of the wave func- 
tions required for the derivation of the bound is intricate. We have calculated them through 
term-by-term integration of the power series expansion of the zero-energy wave functions and 
additionally determining the integration constants that are not generated by this process. 

We define the causal range as the minimum value of the interaction range consistent with 
the causality bound. In effective field theories with contact interactions such as the halo 
EFT, the natural momentum cutoff is of the order of the inverse of the causal range. If 
the natural cutoff is not known from other considerations, its size can be estimated from 
the causal range. If the momentum cutoff used in a calculation is too high, then problems 
with convergence of higher order corrections can appear. Our results can thus be viewed 
as a guide for improving the convergence of EFT calculations with contact interactions. In 
lattice simulations of halo EFT, the lattice spacing should not be taken smaller than the 
causal range. 

We have analyzed the causal ranges for a variety of systems ranging from proton-proton 
scattering to alpha-alpha scattering. Our results for causal ranges in different partial waves 
in these systems typically vary by factors of 2-3. The precise values are quite sensitive 
to small uncertainties in the effective range parameters. In channels with a large negative 
effective range the causal range is very close to zero, which implies that causality provides 
almost no constraints on the range of the interaction in this case. Thus the causal range 
provides a good order of magnitude estimate of the range of interaction, but drawing more 
quantitative conclusions about the structure of the underlying potentials is difficult. 

After an analytic continuation to the bound state regime, the integral relations for the 
causality bound can also be used to derive a model-independent expression for the ANC for 
shallow bound states. If the state is a two-body halo state, the relation can be used to extract 
the ANC from low-energy scattering parameters. Up to higher order shape parameters, our 
relation is equivalent to the one previously derived by Sparenberg et al. in Ref. [22] (see 
Appendix D for details). One difference is that we express the ANC in terms of the binding 
momentum and the effective range rather than the binding momentum and the scattering 
length. We find this form more suitable for the extraction of ANCs from scattering data 
since the effective range is typically more precisely determined than the scattering length 
for shallow states. Moreover, extracting the effective range in a self-consistent fit from the 
scattering data that reproduces the correct separation energy improves the stability of the 
extraction. 

Finally, we illustrate our relation by extracting the ANCs of the excited 2"*" and 1~ states 
in ^^O from a— ^^C scattering data. Compared to previous extractions [HI HH |19] , our values 
are generally smaller. Whether this difference is physically significant requires further study. 
The application of our relation to other shallow cluster states and a benchmark against model 
calculations would also be very interesting. 
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Appendix A: Calculating the constant terms in the causality bound function 



To determine the constants Zi introduced in Sec. Ill D we consider the explicit form of 
g{p,r). From the results of BoUe and Gesztesy [SU], we have'' 

gip,r) =Afeip) ■ 7log(|7|r) ■ f{p,r) 

{ oo 



n=0 



2£ 



+ Re — ^e-'^^-5^rf,,„(p)r" , (Al) 



ra=0 



where 



{2p) 
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r(2£ + 2)2 11 

s=l 



(A2) 



ai,n{p) 



-r(2/ + 2) 



r(n + l)r(n + 2£ + 2) 

r(2/ + 2) 



(2ip)" Wis + l-iri)- [^(n + 1) + iP{n + 21 + 2)] , (A3a) 



s=l 



T{n + l)T{n + 2e + 2) 



(2iprllis + i-iv)-J2-^ 



and 



1 



T(n + 



yy(2ip)"n 

' s=l 



S — 



1 — it] 



s-21-1 



(A3b) 



(A4) 



With this result and the appropriate expression for (70 (r) from Eq. (53) or (54), one can 
use the following procedure to calculate the Z^, z.e., the terms of order in the Wronskian 



1. Note that 



W[g,g^]{r)=p'W[g2.g,]{r) + 0{p') 



(A5) 



and calculate this Wronskian using a truncated version (including terms up to the 
order 2£ + 1 in r is sufficient) of g{p^ r) as given in Eq. ( Al ). 



2. From the result, extract the terms that are of the order r°. 



''One gets this form from Eq. (B6I in Appendix [B] by inserting rt = 2£ + 3, 77 = ^/{2p) and assuming that the 
momentum p is real and positive. 
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3. From that expression then extract the terms that are of the order p^. They constitute 
the (9(r°)-contributions in a series expansion of W[g2, fi'oK'") which cannot be obtained 
from a term-by-term integration of goir)"^- 

With a computer algebra software, this prescription is straightforward to implement. The 
, 2 are shown in Table [T] in Sec. HID 



results for 



0, 



At this point we remark that there is a recent publication by Seaton [51j discussing 
Coulomb wave functions that are explicitly analytic in the energy (or the momentum 
squared). In principle, it is possible to use these results to get explicit expressions for 
/2(r) and (72 (r) (in addition to the already known zero-energy functions), and then simply 
calculate all the Wronskians directly. However, the analytic irregular Coulomb function de- 
fined by Seaton is slightly different from our g{p,r). More importantly, not all coefficients 
needed for the expansions in energy are given explicitly. Finally, writing everything in terms 
of the wave functions of Bolle and Gesztesy paves the way for a generalization of the results 
presented in the current paper to an arbitrary number of spatial dimensions. 

Knowing now what the constant terms Zi that come from W[g2, go]{r) should be, it is 
possible to write down explicit expressions for the causality bound functions bf'{r). To do 
that we use that the antiderivative of the right hand side of Eq. (52) can be expressed 
in terms of (generalized) hypergeometric functions pFg and Meijer G-functions G^^. The 
only additional point to be taken into account is that the antiderivative of go{rY in general 
includes a constant term that is different from the desired Z^. Hence, one has to determine 
this term and add another constant such that the sum is exactly equal to the Zi given in 
Table |l| For the antiderivative of goirY that we are using in the following, up to a minus 
sign the correction term is just the AZ^ introduced in Sec. |vj^ 



a. Repulsive case, 7 > 



get 



The i = result for a repulsive Coulomb potential is given by Eq. (63). For i = 1, we 
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(A6) 



G^]^ 1 47r 



-3,-2,0,3; 54' 



^In other words, our choice for the antiderivative corresponds to using dr'gQ{r')^, at least for the case of 
a repulsive Coulomb interaction, where this integral is convergent. 
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and for ^ 



2 the result is 



687' 



77V^ 1F2 ( ^; 5, 10; 4jr] - 1087^^ .F^ 4, 9; 47r 



+I2967V" 1F2 (^-; 3, 8; 47rJ - 120967r 1F2 (^-; 2, 7; 47r 
+ 108864 1F2 Q;l,6;47r') - 108864 2F3 Q, 2; 1, 3, 6; 47r 



4^2 



K)"'G|n47r 



0,5,-5,-2 



0r7^r^ 
3600 



( 47r 



-5,-2,0,5 



+ 



7 



21600 



b. Attractive case, 7 < 



For an attractive Coulomb potential, the 
result reads 
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result is given by Eq. (64). The 
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and for ^ 



2 one finds 
50r2 
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Appendix B: The Coulomb wave functions of Bolle and Gesztesy 

BoUe and Gesztesy define the Coulomb wave functions 

r) = r^+™e-'P"iFi (| + m - 1 + 2m; z) 



and 



(P, r) = ^pt^L,''^ (2ipf"^ ri+™e-'^^f/(^ + m - 1 + 2m; ^) , 
1 (2m + 1) ' 



(Bla) 



(Bib) 



where, in our case, n = 2£ + 3 and m, k, z are as defined in Eq. (11). Using the formula^ 

2V^'?/2 \Y{i+l + ir])\ 



a 



ri,i 



{2i+l) 



(B2) 



Eq. (5.1) in Ref. [26], one finds that 
and 

Gi°)(p,r)=/C,,, 



(p)/ 



(B3a) 



(B3b) 



It is shown in Ref. [30] that Fn\p-, r) is analytic in p'^. Furthermore, from Eqs. (3.16), (3.17) 
and (4.1) in that paper it follows that 



(B4) 



^From Eq. (3.1) in Ref. ^25] one directly sees that this is consistent with Eq. (27) in Sec. II B 
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where he{p) is the function defined in Eq. (26), and where Gn\p,r) is analytic in p^. The 
functions f{p,r) and g{p,r) that we introduce in Sec. IIIB are exactly the analytic wave 
functions defined above, i.e., 

f{p, r) = Fi'\p, r) and ~g{p, r) = (P, r) . (B5) 

It should be noted that Lambert [32] already defines Coulomb wave functions that are 
analytic in p^, but only for the most common case of three spatial dimensions. BoUe and 
Gesztesy extend Lambert's results to an arbitrary number of dimensions > 2. 

More importantly, at least for our application, in Eq. (4.3) of Ref. [30], BoUe and Gesztesy 
give an explicit expression for Gn\p,r). Since there are two typos in their original equa- 
tion, we quote the whole expression for completeness. Slightly altering the notation to 
match our conventions, we have 
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(B6) 



valid for any odd n > 3 and where the sum over s in the fourth line is defined to give zero 
for n = 3 and A; = 0. 



Appendix C: Numerical calculations 

In order to check our relations and to get a better understanding of the values for the 
causal range, we consider some explicit examples numerically. 

By cutting off the singular parts of the potential {i.e., the Coulomb potential and the 
angular momentum term for £ > 0) at very small distances, it is a simple task to numerically 
solve the radial Schrodinger equation Q in configuration space. From the radial wave 



'The \^/k in the first fine should be i7/(2fc) and the {q + 1) in the last line should be V{q + 1). 
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functions one can extract the Coulomb-modified phase shifts by looking for a zero at some 
large {i.e., much larger than the range R of the short-range potential) distance, 



w 



ip) 



(ro) =0 , ro>i?, 



and then calculating 



cot6e{p) 



Ft{p){ro) 



For the simplest case of a local step potential, 



V{r, r') = V3tcp(r) ■ 6{r - r') = Vq e{R - r) ■ 6{r - r') , 



(CI) 



(C2) 



(C3) 



one can of course also obtain the phase shift directly by matching the wave functions at 
r = R. The effective range parameters are then obtained by repeating the calculation for 



several (small) momenta and fitting Eq. (21) to the results. 

In order to test Eq. (44) directly one needs the wave function to calculate the integral 



IT 

I dr' 




2 r 
= lim / dr' 


(p)/ '\ 
wf [r ) 


Jo 









Even if we do not actually take the limit p — )■ but rather just insert some small po 
(in units of an arbitrary inverse length scale), we find that the relation 



bf{R) 



R 



dr' 



(C4) 
= 0.1 

(C5) 



is typically fulfilled to better than one percent accuracy for the simple step potential defined 
above. 

For illustration, in the following we choose units where the radial distance is measured in 
fm. The potential range is set to 1 fm and its strength is measured in MeV.^^ Furthermore, 
the reduced mass and Coulomb parameter are set to the values for the proton-proton system, 
z.e., 2fi = mN^ 940 MeV and 7 = 7p_p ^ 0.035 fm'^ 

In Figs. [3] and [4] we show the results (for £ = 0,1,2) for both repulsive and attractive 
step potentials. Quite interestingly, the £ = causal range for the repulsive potential stays 
at zero (meaning that one could reproduce the underlying values of the scattering length 
and the effective range even with a contact interaction) until a potential strength of about 
100 MeV. For higher partial waves the causal range takes a nonzero value for much weaker 
potentials, but the rise is less steep. In general, it is remarkable that the causal range is 
typically considerably smaller than the actual potential range {R = 1 fm). 

For attractive potentials the causal range grows much faster as the potential strength 
(now negative) increases. In contrast to what one might expect, no special features are seen 
in the causal ranges as the potential becomes strong enough to support a new bound state 
close to threshold, i.e., when there is a pole in the scattering length parameter. 

To conclude this section, we show the general dependence of the causal range on both 
the scattering length and the effective range, which has the advantage of not depending on 



^^Note that with these conventions, the quantity that is used in the numerical calculation is vq = 2^Vq/ {hc)^ , 
where he w 197.33 MeV • fni is used for the unit conversion. 
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FIG. 3: Causal range for a repulsive step potential and 7 = 7p-p. 
0.8 




-150 -100 
Vq (MeV) 



FIG. 4: Causal range (thick lines) for an attractive step potential and 7 = 7p-p- The thin lines 
show the corresponding scattering lengths. 



a certain model potential. For illustration, we again measure distances in fm and set the 
Coulomb parameter to the value of the proton-proton system. In Figs. |5] and |6] we show 
the results for i = and i = 1. For negative effective range r'^ , the causal range stays 
essentially zero. For positive effective range, it increases as the absolute value of a*" becomes 
larger. If one gradually turns off the Coulomb interaction by letting 7 — )■ 0, the i = 1 plot 
stays almost unchanged, whereas the i = result remains qualitatively the same, but with 
a much steeper rise in the quadrant where a*" > and r'-^ > 0. 
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£ = 0, 7 = 7p-p 




Rc (fm) 




FIG. 6: Causal range for 7 = Jp-p and £ = 1 in dependence of and , measured in fm^ and 
fm~^, respectively. 

Appendix D: Equivalence to relation by Sparenberg et al. 

In this appendix, we demonstrate the equivalence of our ANC relation to the one derived 
by Sparenberg et al. [22]. 
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Their ANC relation, translated to our notational convention, reads 



1^1 



K 



£+1 



(Dl) 



The tilde on the is there to indicate that Sparenberg et al. use a convention for the 
Coulomb-modified effective range expansion that is slightly different from ours (apparently, it 
comes from Ref.[52]). Combining equations in Ref. [22] and again matching to our notation, 
one finds that 



£!2 



n(^' + r,') \CI,V''^' cot Uv) +7P'' Hp) 



s=l 



(D2) 



For S- waves this is equivalent to our Eq. (21), but in general the left-hand sides of the 



expansions differ by an overall factor [2^i\/{2i + 1)!]^. More precisely, 



a? 



(2£ + l)! 



(2£+l) 
i\ 2^ 



r? , etc. 



(D3) 



The final step of the derivation by Sparenberg et al. eliminates the effective range in 
favor of the scattering length. Without invoking this final step their relation reads 



\A\ 



2^(2£ + 1)! 
r{i+l + fi) 



+ 



7 



2£-l 



3(2£+ 1)!2 



^1/2 



(D4) 



where we are now using our convention for the effective range exp ansion. With the definition 
of Crj/ from Eq. (76) and the values for AZi 

i = o[ 



from Table 



III 



one sees that at least for 

, 2 this is exactly equivalent to our Eq. (84) with the ^(k^) set to zero. In order to 
prove the equivalence for arbitrary i one would need a general expression for AZ^. This, in 
turn, requires knowledge of the constant terms in W[g2, go]{f) for arbitrary i. It would thus 
probably be more interesting to turn the argument around, i.e., take the equivalence for 
granted and derive from it a general expression for the constant terms in the Wronskians. 
The only additional ingredient one would need for this procedure is a general series expansion 
for the Meijer G-functions that arise from the integral of goir)"^- 
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